Statistical Physics of Elasto-Plastic Steady States in Amorphous Sohds: Finite 

Temperatures and Strain Rates 



O 

1—5 

oo 



o 

s 



I 

o 



> 

m 
i> 
cn 

o 
o 



X 



Smarajit Karmakar, Edan Lerner, Itamar Procaccia and Jacques Zylberg 
Department of Chemical Physics, The Weizmann Institute of Science, Rehovot 76100, Israel 

(Dated: June 21, 2010) 

The effect of finite temperature T and finite strain rate 7 on the statistical physics of plastic 
deformations in amorphous solids made of particles is investigated. We recognize three regimes 
of temperature where the statistics are quahtatively different. In the first regime the temperature is 
very low, T < Tctoss{N), and the strain is quasi-static. In this regime the elasto- plastic steady state 
exhibits highly correlated plastic events whose statistics are characterized by anomalous exponents. 
In the second regime Tcross{N) <T < rmax(7) the system-size dependence of the stress fluctuations 
becomes normal, but the variance depends on the strain rate. The physical mechanism of the cross- 
over is different for increasing temperature and increasing strain rate, since the plastic events are 
still dominated by the mechanical instabilities (seen as an eigenvalue of the Hessian matrix going 
to zero), and the effect of temperature is only to facilitate the transition. A third regime occurs 
above the second cross-over temperature rmax(7) where stress fluctuations become dominated by 
thermal noise. Throughout the paper we demonstrate that scaling concepts are highly relevant for 
the problem at hand, and flnally we present a scaling theory that is able to collapse the data for all 
the values of temperatures and strain rates, providing us with a high degree of predictability. 



I. INTRODUCTION 

The simulational investigation of elasto-plastic steady 
states in amorphous sohds concentrated in recent years 
on athermal and quasi-static conditions (AQS) [1-5]. 
The reasons for this are manifold: first, in athermal con- 
ditions the separation between elastic regimes and plas- 
tic events is clear cut. Second, the mechanism for plastic 
events is purely mechanical, and it can be understood 
entirely on the basis of the underlying inter-particle po- 
tentials and dynamics [6, 7]. Third, the simulations ex- 
hibit the existence of fascinating, highly correlated plas- 
tic events whose statistics abound with anomalous ex- 
ponents that are not fully understood. On the other 
hand, current models of elasto-plasticity tend to propose 
mean-field type theories in which plastic events are un- 
correlated and are not linked to mechanical instabilities, 
but are rather assumed to be triggered by purely thermal 
fluctuations [8-13]. If so, the interesting findings in the 
AQS regime might be entirely irrelevant for "normal plas- 
ticity" at higher temperatures and strain rates. The aim 
of this paper is to examine this issue with some care. To 
this purpose we will explore here the effect of tempera- 
ture and of strain rates, separately and together, to assess 
whether indeed the existence of finite temperatures and 
strain rates obliterate the relevance of the rich plethora 
of findings at the AQS conditions. Our conclusion is that 
this is far from being true [14, 15]. In fact, even when 
temperature or strain rates are enough to make the stress 
and energy fluctuations 'normal' (in a sense made precise 
below) the role of mechanical instabilities is still crucial. 
We thus need to reveal the mechanisms responsible for 
the cross-over in statistics and to propose a theoretical 
framework for the discussion of the statistical physics of 
elasto-plastic steady state taking into account all the es- 
sential ingredients. When done, as described below in 
this paper, we will own a description that unifies the 



statistical physics below and above the crossovers from 
highly correlated to uncorrelated stress and energy fiuc- 
tuations. 

The structure of this paper is as follows. In Sect. II we 
present the model glass former and the details of the nu- 
merical procedures employed throughout this study. In 
Sect. Ill we present a short review of the essential re- 
sults pertaining to elasto-plasticity in amorphous solids 
in AQS conditions. In Sect. IV we consider the effect 
of temperature on elasto-plasticity with very slow strain 
rates (where very slow means rates below the cross-over 
to independent stress and energy fluctuations due to high 
strain rates). We show that the plastic events are as- 
sisted by thermal fluctuations but they are nevertheless 
dominated by mechanical instabilities. The cross-overs 
from highly correlated plastic events to independent fluc- 
tuations is considered in Sect. V, where the different 
mechanisms for these cross-overs due to thermal fluctua- 
tions and due to high strain rates are treated separately. 
Throughout the paper we use scaling concepts to orga- 
nize, compactify and exhibit the physics in its neatest 
form. This approach culminates in Sect. VI where we 
present the scaling function theory that unifies the statis- 
tical physics in the elasto-plastic steady state before and 
after the thermal-dominated and strain-rate-dominated 
cross overs. Sect. VII offers a summary of the paper and 
concluding remarks. 



II. MODEL AND NUMERICAL METHODS 

A. System Definitions 

Below we employ a model glass-forming system with 
point particles of two 'sizes' but of equal mass m in two 
and three dimensions (2D and 3D respectively), interact- 
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ing via a pairwise potential of the form 



2i 



(1) 

where rjj is the distance between particle i and j, £ is the 
energy scale, and Xc is the dimensionless length for which 
the potential will vanish continuously up to q derivatives. 
The interaction Icngthscalo Ay between any two particles 
i and j is = l.OA, A^ = 1.18A and A^ = 1.4A for two 
'small' particles, one 'large' and one 'small' particle and 
two 'large' particle respectively. The coefficients C2^ are 
given by 



(-1) 



(fc + 2g)!! 



(2g-2£)!!(2£)!! (fc - 2)!!(fc + 2^) 



-(fe+2£) 



(2) 



We chose the parameters Xc = 1.385, A: = 10 and 
q = 2. The unit of length A is set to be the interac- 
tion length scale of two small particles, and e is the unit 
of energy. Accordingly, the time is measured in units 
of = \Jra}? je. The density for all systems is set 
to be A^/F = 0.85A^^. The glass transition tempera- 
ture is Tg K. 0.46£/A;b, defined here by the condition 
7"a(T'g) = lO^r*, where Tq, is the structural relaxation 
time. 



B. Methods 

The work presented here is based on three types of 
simulational methods. The first type corresponds to the 
athermal quasi-static (AQS) limit T — )• and 7 — )• 0, 
where 7 is the strain rate. AQS simulations had been 
extensively analyzed recently [1-3, 5, 16, 17] as a tool for 
investigating plasticity in amorphous systems. In AQS 
simulations one starts from a completely quenched con- 
figuration of the system, and applies an afhne simple 
shear transformation to each particle i in our shear cell, 
according to 



rix + TiyS-y 



(3) 



in addition to imposing Lees-Edwards boundary condi- 
tions [18], and (57 = 7 — 70 is a small strain increment 
from some reference strain 70. The strain increment Sj 
plays a role analogous to the integration step in standard 
MD simulations. We choose the basic strain increment 
step to be (57 = 5 X 10~^ for all system sizes simulated, 
and sample each plastic event using strain increments 
of at most 2 x 10~^ using the backtracking procedure 
described in [5]. The afRne transformation (3) is then 
followed by the minimization of the potential energy un- 
der the constraints imposed by the strain increment and 
the periodic boundary conditions. Wc chose the termina- 
tion threshold of the minimizations to be iVit/l = 10~^, 



for every degree of freedom Xi. Our method for locating 
saddle points is explained in Subsect. IV B. 

The second simulation method employs the so-called 
SLLOD equations of motion [18]. For our constant strain 
rate 2D systems, they read 

hx = Pix/m + jny , 

'^iy ~ Piy/'rn , 
Vix ~ fix "jPiy ? 



Piy — fil 



(4) 



We use a leapfrog integration scheme for the above equa- 
tions, and keep the temperature constant by employing a 
modification of the Berendsen thermostat [18], measuring 
the instantaneous temperature with respect to a homo- 
geneous shear flow. The modification implies that we 
randomly re-partition the system into subsets of about 
500 particles, and utilize a set of Berendsen factors, with 
a different factor for each subset (instead of just one in 
the standard algorithm). This modification was found 
necessary in order to reduce finite-size effects due to sub- 
extensive statistics. The integration time step was cho- 
sen to be 6t = 0.005, and we set the time scale for heat 
extraction at tt = 5. Or*. 

The third simulational method [19] was employed to 
study systems at athermal conditions, but at finite strain 
rates. This method utilizes the SLLOD equations of mo- 
tion (4), with the addition of total momentum conserving 
damping forces, such that the force on the i'th particle 
is given by 



TT '—f. 



where is calculated from the pair potential 
Dij is given as 

D,, = l-2in,/r,)^ + in,/r,f 



(5) 



and 



(6) 



which vanishes smoothly at = Tc, where Tc is the cut- 
off of the pair potential (f)ij . For integrating the equation 

of motion we will use a slightly modified version of the 
Velocity- Verlet algorithm defined below [20] 



V^{t)St■ 



Vi{t + 6t) = Vi{t) + \ {fi{t) + fi{t + 5t)) 5t . 



(7) 



Below we need to determine whether a thermal stressed 
system still resides in an original local minimum of the 
athermal system or whether it had jumped to the basin 
of attraction of another local minimum. This is done 
by taking a given configuration and minimizing its po- 
tential energy until it hits the minimum. For a system 
that is stressed at a finite temperature one observes plas- 
tic events before the mechanical instability threshold is 
reached. By minimizing the energy every 100 times steps 
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in the simulation below we can identify such transitions 
by finding that the original minimum is no longer cap- 
tured and a new one replaced it. 

For completeness we also carried out simulations in 
three dimensions, using the same binary mixture with the 
same interaction potential as in 2D but with a density 
N /V = 0.81. Some result from these simulations are 
discussed below. 



III. REVIEW OF ELASTO-PLASTICITY IN AQS 
CONDITIONS 

We consider amorphous solids in the limit of zero tem- 
perature T — > 0, subjected, say, to shear deformation 
at vanishing low strain rates 7, with 7 our parametriza- 
tion of the imposed deformation, see below. An amor- 
phous solid in the athermal limit must satisfy the follow- 
ing conditions [3]: {i) the notion of solidity requires that 
all the eigenvalues of the Hessian matrix Ti^ = g^.g^. 
are strictly positive, (ii) The amorphous nature of the 
considered systems is guaranteed by demanding that the 
mismatch forces Hi = F are non-zero and uncorre- 



lated, i.e. ^ 0, (SiSj 



(Hi) The limit T 





implies that the system always resides in a local minimum 
of the potential, with the forces 



dU 
dxi 







(8) 



at all times. 

In AQS simulations the potential energy [/ is a function 
of the imposed strain, parameterized by 7, and of the 
particle coordinates £Ci(7), U = ?7({a;i(7)}, 7). Below 
we derive the explicit coordinate dependence on 7; we 
consider deformations via parameterized transformations 
on the particle coordinates J? (7) — X ^ •jh (not to be 
confused with the Hessian H): 



H Xi 



Ui 



(9) 



where the non-affine coordinates Ui are additional dis- 
placements that assure that the zero-forces constraint (8) 
is fulfilled; total derivatives with respect to strain in the 
athermal limit should thus satisfy the zero-forces con- 
straint. They are carried out via 



d d dui d 
d'j dj d'y dui 



d dui d 
97 c?7 dxi 



(10) 



where the second equality results from from Eq. (9), and 
here and below repeated indices are summed over. The 
evolution of the non-affine coordinates Ui can be explic- 
itly derived by requiring that ^ = 0: 

dfi _ dfi dfi duj 



d'j dj dxj dj 



= . 



(11) 



We refer to the full derivatives of the relaxation coordi- 
nates with respect to strain as the non-affine velocities 



Vi = Inserting the definitions of the Hessian Hij 

and the mismatch forces in Eq. (11), we obtain 



Sii -\- Hij ' Vj — . 



(12) 



From here, the equation for the non-affine velocities {vi} 
is obtained by inverting (12): 



(13) 



With an equation for the non-affine velocities Vi, full 
derivatives with respect to strain can be written as 



d _ d 
d'j d"f 



d 
dxi 



(14) 



and the equation of motion for the coordinates Xi(-^) be- 
comes explicitly available: 



dxi 

— — = h ■ Xi + Vi . 

07 



(15) 
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FIG. 1: A typical stress vs. strain curve in a system of 4096 
particles under simple shear deformation in two dimensions 
obtained in the athermal quasistatic limit. Every elastic (re- 
versible) increase in stress is followed by a sudden plastic (ir- 
reversible) drop in stress. 

The above analysis of the AQS dynamics is valid as 
long as all the eigenvalues of the Hessian are positive; this 
condition breaks down eventually upon increasing the ex- 
ternal strain, when reversible elastic branches are ter- 
minated by mechanical instabilities, as demonstrated in 
Fig. 1. These mechanical instabilities are associated with 
the vanishing of an eigenvalue of the Hessian, when the 
local minimum at which the system resided develops an 
unstable direction in coordinate space, which the system 
follows while undergoing an abrupt drop in energy. The 
system then descends down the potential energy land- 
scape until finding some other mechanicaly stable local 
minimum. Associated with this drop in energy, denoted 
in the following by AU — J/initiai ~ %nai, is a change in 
stress, denoted in the following by Act = CTinitiai — o-Rmd, 
the change in stress is not strictly constrained to be posi- 
tive - in small systems one may encounter mechanical in- 
stabilities that result in an increase rather than a decrease 



4 



of stress. In this section wc will provide a brief review 
of the physics of deformed systems approaching mechan- 
ical instabilities, followed by a review of the statistics of 
the flow events in the steady flow state. We will demon- 
strate in the next section that the mechanical instabilities 
remain highly relevant also in the second temperature 
regime Tljross 



A. Mechanical instabilities 

As mentioned above, upon increasing the external 
strain, an eigenvalue of the Hessian eventually vanishes 
at the inevitable onset of a mechanical instability. We 
denote the vanishing eigenvalue as Ap, its corresponding 
eigenfunction as and the strain value at which the 



upcoming instability occurs as 7p, such that 
Ap — >■ as 7 — >■ 7p . 



(16) 



Returning to Eq. (13), we write the non-affine velocities 
in a normal mode decomposition [1], 



/ (fe) — 



k 



where Afe is an eigenvalue of the Hessian and i/'j'''' its 

corresponding eigenfunction: "Hij ■ ipj = Xk'^j . 
Ap — > 0, the above sum will be dominated by the diverg- 
ing term, i.e. 



As 



Af 



as 7 — >■ 7p 



(18) 



We now calculate derivatives of the potential energy with 
respect to strain in the vicinity of a mechanical instabil- 
ity, say, at 7p; when possible, wc will only keep terms 
that consist of contractions with the non-affine velocities 
Vi, since these diverge as the system is strained towards 
7p, cf. Eq. (16) and Eq. (18), hence they will domi- 
nate over any regular terms. With V = L'^ denoting the 
volume of the system with linear size L, the stress a is 
deflned via 



^^_dU _ dU 



dU 
dxj 



where the second equality stems from the constraint of 
zero forces (8). Since the non-afflne velocities do not 
appear in the above expression for the stress, it is always 
regular, even as 7 — >■ 7p. Taking another total derivative 
of the potential energy with respect to strain, keeping 
only highest order terms in : 



dxjdj 



Vi • a,- ~ - 



Ap 



(20) 



with a = {ipj^^ ■ 3j)^. Taking a last total derivative of 
the potential energy with respect to strain requires an 



expression for with the notation Tijk = g^S^idxi ' 
an expression can be obtained by writing 



dj"^ dj 



— {Si + nij-vj) 



= + Tijk ■ VjVk 



(21) 



and inverting for Keeping the highest order terms 
in we find that in the vicinity of 7p, 



dvj 



-'Hi/ ■ {Tjke ■■ vkve) , 



thus 



d^U ^ dVj ^ 
d'y^ dj 



(22) 



(23) 



with b = {^P-Stf Tijk-tpPil^f^tpT^- Combining re- 
(17) suits (20) and (23), we arrive at the differential equation 



d_ 

dj \Xp 



(24) 



The solution of the above difFcrctial equation, together 
with the boundary condition Ap(7p) =0, is [6], 



Ap ~ V 7p — 7 



(25) 



We should emphasize at this point that the potential 
energy should, in principle, be written as a series ex- 
pansion in terms of components of the strain tensor 
e = ^{H'^ H I) [7]. Then, additional terms may ap- 
pear in potential energy derivatives, given a parametriza- 
tion of H{j). Here, for the sake of simplicity, we di- 
rectly utilize a parameterized notion of deformation via 
the parameter 7. Since we were only interested in sin- 
gular terms near mechanical instabilities, our results are 
equally valid. 



B. Flow events statistics 

When subjecting an amorphous solids to external shear 
strain, it tends to set up an elasto-plastic steady flow 
state after a short transient of a few percent strain. In 
the steady flow state, the statistics of the energy drops 
AU, the stress drops Act and the strain intervals between 
successive flow events A7 become stationary. In partic- 
ular, one finds that the averages of these quantities obey 
the following scaling relations 



(AC/) 
(A7)' 



(Act) ~ sN^ 



(26) 
(27) 
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In Fig. 2 the mean energy drop (A/7) and mean strain in- 
terval (A7) for our model system are displayed, together 
with the scaling laws (26). In the upper panels we show 
results in two dimensions and in the lower panel in three 
dimensions, In the present model we find that a « 1/3 
and 13 « -2/3 in both 2D and 3D. 




FIG. 2: Mean energy drop (At/) (left panels) and mean strain 
interval (A7) (right panels) as functions of system size, mea- 
sured in AQS simulations of steady plastic flow of the system 
described in Sect. II. Upper panels: two dimensions. Lower 
panels: three dimensions. The continuous lines represent the 
scaling laws (26) and (27). The scaling exponents are the 
same in 2D and 3D. 

A scaling relation a — /3 — 1 follows from the average 
energy balance equation, cf. [5] 

aY{Aa)V/fi = (AU) , (28) 

where cry- is the yield stress (the mean stress in the steady 
state flow state, in the AQS limit), and fi is the shear 
modulus. 



C. The scaling of the variance of stress fluctuations 

For our purpose of distinguishing clearly between dif- 
ferent regimes of elasto-plastic statistical physics it is ad- 
vantageous to measure the properties of stress fluctua- 
tions, and in particular of the variance. This quantity 
will reflect the change in physics in the different regimes. 
In AQS conditions the variance of stress fluctuations 
{Sa"^) = {{a — (JyY) also exhibits anomalous scaling laws 
with the system size. In particular we find 

(5^^2) _ ^ ^29) 

with « —0.4 in the studied model, cf. Fig 3. One 
should notice the difference between the exponent char- 



acterizing the N dependence of \J {Sa'^) and of the ather- 
mal mean stress drop (Ac), in the sense that 9 f3. This 
difference is due to very strong correlations between elas- 
tic increases and plastic drops. 
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FIG. 3: Stress fluctuations (Sa^) of the steady flow state of 
the studied model. The continuous line represents the scaling 
law (29). 



IV. THE MECHANISM OF THERMALLY 
ACTIVATED PLASTICITY 

As explained above, in our thinking about thermal 
effects on elasto-plasticity in finite systems strained at 
finite strain rates and finite temperatures there exist 
three temperature regimes in which the dynamics and 
the statistics of plastic events are qualitatively differ- 
ent. In this section we investigate the second regime 
TcrossiN) <T< Tinaxli) whcrc Tcross(A^) IS systcm size 
dependent as explained in the following subsection. 

A. The first crossover in statistics 

The flrst cross-over in statistics from avalanches to 
independent statistics occurs at the temperatures T « 
Tcross where the cross-over temperature is estimated as 
follows [15]: during a plastic drop the energy released 
spreads out quickly in the system on the time scale of 
elastic waves. Thus every particle shares an energy of 
the order of eN"/N = eN^. On the other hand the typ- 
ical scale of thermal energy per particle is T (in units of 
Boltzmann's constant). We thus expect thermal effects 
to start overwhelming the statistics of athermal plastic 
events when 

^ Tcross . (30) 

This equality will hold when the system size L = ^2, 
where (^2/^)"* = N. Substituting the last equality in Eq. 
(30) and then solving for ^2 we find 

6/A = [nross/lf^ ■ (31) 

Note that Eq. (30) implies (since /? < 0) that the change 
in statistics occurs at a temperature that decreases when 
the system size increases, meanings that AQS statistics 
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FIG. 4: The mean energy drop in plastic events as a func- 
tion of system size for different temperatures, see inset. The 
straight line represents the AQS scaling of mean energy drops. 



will pertain only for small systems. We also understand 
the meaning of the cross-over due to thermal effects: the 
thermal agitation reaches just the necessary level to com- 
pete with the stored elastic energy per particle. The 
avalanches are made of the primary plastic instability 
triggering all the other, close to instability regions, to 
flip in tandem and relax. These ripe regions which are 
sufficiently close to instability to react to the primary 
instability, are all destroyed by the thermal fluctuations 
such that the primary instability remains naked, turn- 
ing the statistics of the plastic events from anomalous to 
normal. 

This picture can be demonstrated directly by measur- 
ing the energy drop in a plastic event in quasi-static con- 
ditions but at different temperatures. The measurement 
is done by stopping the thermal molecular dynamics sim- 
ulation, followed by quenching the system to a very low 
temperature of T = 10^"^ on a time scale of lOOr^. This 
procedure allows for the completion of any plastic ac- 
tivity. Then the potential energy is minimized and the 
AQS scheme is employed to measure the energy drop in 
the flrst upcoming mechanical instability. This method 
allows us to probe directly the consequences of the ther- 
mal agitation on the ability of the system to undergo an 
avalanche. The results of such measurements are shown 
in Fig. 4. We see that the potential AQS energy drops are 
sapped out by the thermal agitation. Only at the lowest 
temperature of T = 0.01 the energy drops approach the 
AQS limit for small systems, but even for this low tem- 
perature larger systems cannot come close to the AQS 
limit. This is all in accordance with the estimates in 
Eqs. (30) and (31). We stress that this cross-over occurs 
at quasistatic conditions (but not only) and has nothing 
to do with 7. 

An additional direct demonstration of the first thermal 
cross-over is obtained by measuring the stress fluctua- 
tions. Recall that at AQS conditions these fluctuation 
exhibit anomalous scaling , cf. Fig. 3. At higher tem- 
peratures the data in Fig. 5 indicate a clear cross-over to 
independent stress fluctuations in which {Sa'^) ~ iV~^. 




10'' 10^ 10" 10^^ 10"' L 10' 



FIG. 5: Left panel: the variance of the stress fluctuations 
(multiplied by A'^) as a function of the system size A*' for a 
2D system and for various temperatures, at a strain rate 7 = 
2 X IO"*". Right panel: The scaling function g{x), cf. Eq. 
(36) for the data in the left panel. Note the cross-over for x 
of the order of unity as predicted by Eq. (34). The power 
law decrease at low values of x are in agreement with the 
prediction of « 0.33. The two black hues represent the 
theoretical prediction for the scaling function g{x) for a; ^ 1 
and for a; ^ 1. Note that the full scaling function depends on 
7, and the present one is an approximate version for 7 — >■ 0. 

To capture the temperature and size dependence of the 
variance, and to demonstrate unequivocally the thermal 
cross-over, we first need to separate the thermal (vibra- 
tional) contribution from the mechanical contributions to 
{5cr^). We write 

(Sa^) = {6a^)T + (S^) , (32) 

where {Sa'^)T denotes the thermal contribution which can 
be read from Eq. (10) of Rcf. [22], i.e. 

{6c7^)t « AiT/y . (33) 

For the mechanical part we introduce a scaling function 
which exhibits the thermal cross-over. In other words, we 
propose a scaling function 52(2;) to describe the system- 
size and temperature dependence of the mechanical part 
of the variance. To this aim we introduce the tempera- 
ture dependent natural scale ^2(^^)1 

^2{T) = X[T/e]'/''^ . (34) 

The scaling function 32 (x) is a function of the dimension- 
less ratio L/£^2'- 

{d^){N,T) = s^N^'g^iL/UT)) ■ (35) 

The dimensionless scaling function 172(2^) must satisfy 

32(2;) go; for a; ^- , 

g2ix) — gox'' for X —7- 00 . (36) 

The first of these requirements means that the fiuctuation 
are in accordance with the athermal limit. The second 
means that after the cross-over the fluctuations of the 
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stress become intensive, requiring ( = —d{l + 29). We 
compute C « —0.4 in 2D. Note that in Ref. [15] the same 
scahng function was written in terms of 1/x instead of x. 

We present tests of the scahng function in Fig. 5. Ex- 
amining the right panel of Fig. 5 we see that the ther- 
mal cross-over is demonstrated very well where expected, 
i.e. at values of x of the order of unity. The asymp- 
totic behavior of the scaling functions agrees satisfacto- 
rily with the theoretical prediction which are indicated 
by the black lines. 



It was shown in Ref. [23] that the energy barrier AE{'~f) 
obeys the scaling law 



A£;(7) c(7p - 7)^ 



(41) 
Indeed, 



over very large strain intervals of up to 10 
solving ^ = in Eq. (38) results in the estimation of 

the energy barrier of AE^j) ~ Plugging Eq. (25) 

into this relation, assuming (as is the case in a saddle 
node bifurcation) that 7p is not singular near 7p, we 
obtain 



B. Strain dependent energy barriers 

The aim of this section is to argue that for tempera- 
tures that are not too high the plastic events under exter- 
nal strain are dominated by mechanical instabilities and 
are only assisted by thermal fluctuations. It is crucial 
at this point to clarify what we mean by 'not too high'. 
We will argue that in the elasto-plastic steady state there 
exists a 'typical barrier for thermal activation'. Denoting 
this typical barrier by AGtyp we estimate Tmax by com- 
paring the typical escape time over the typical barrier 
r = exp (AGtyp/T) to (Ty/(/x7). The ratio T^j/ay 
measures the effect of the strain rate and when it ex- 
ceeds some threshold, T^ax should become a decreasing 
function of 7. We estimate AGtyp below and address 
the effect of strain rate in more detail in the next two 
Sections. The present considerations apply for the tem- 
perature range Tcross <T< Tmax- 

Consider an athermal amorphous solid close to a me- 
chanical instability at jp; there, an eigenvalue Xp of the 
Hessian vanishes, as seen in Eq. (25), Xp ^ -\/7P — 7. 
Denoting ■j/;^^-' the eigenvector associated with the van- 
ishing eigenvalue A p , we define the reaction coordinate s 
as the displacement of the system from the minimum of 
U at a;'" in the direction of ip^'^\ i.e. 



(37) 



We expand the potential up to third order in s: 

U{s) ~ [/(a^") + iApS2 + iTpS^ + 0(S4) , (38) 



where 7p 



Vr^r^- The re- 



lation (38) is valid for 7 — 7p; this form of the poten- 
tial energy results in the existence of a saddle point at 
s = —^f^, which translates to the real-space position of 
the saddle point at 
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A£;(7) - - (7. 



7)^ 



(42) 



Note that the scaling Ap ^ y^jp — 7 holds only ex- 
tremely close to 7p (only vanishingly close to 7p as — 
00), while the scaling of the energy barrier, Eq. (41), 
holds along strain scales that are larger by orders of mag- 
nitude, [23]. We show below that the range of strain on 
which the scaling law (41) holds, as well as the pre-factor 
c in (41) are independent of system size. 

To this aim we choose another reaction coordinate de- 
noted s as the displacement of the coordinates x from 
the minimum at a;™, but this time directed towards the 
saddle point, i.e. 



Xi Xa 



sz. 



where 



Ix" — X' 



(43) 



(44) 



Since it is not possible to know the direction 2(7) far 



---[/ 

^min 
^saddle 



(39) 



We define the energy barrier Ai?(7) as the difference be- 
tween the potential at the saddle point /7(a;*), and the 
potential at the minimum [/(a;™), i.e. 



A£;(7) = U{x') ~ C/(a;") 



(40) 



FIG. 6: Illustration of the procedure of finding the potential 
barriers. Since the exact direction of the barrier z is unknown, 
the initial guess brings the system to a state in which U > 
twaddle. The point s* is the displacement along the estimated 
z at which the system leaves the original basin of attraction. 
At s* the saddle point x" can be detected by minimizing the 
square of the gradient function jViC/j^. 
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away from 7p, we begin our measurements of Ai?, for 
various system sizes, very close to 7p, where 2(7 — >■ 
7p) — ?► ip^^K and set z — ■0^^-' as the initial guess for 
z. Then, to find the energy barrier, we displace the sys- 
tem along z by small increments of the reaction coor- 
dinate s, and minimize the potential energy after each 
displacement. For small displacements the minimization 
brings the system back to a;™; however, at some displace- 
ment s*, the system does not return to the minimum at 
a;™, but rather finds a different minimum, see Fig. 6. At 
the displacement value at which the system leaves the 
original basin of attraction, we minimize the square of 
the gradient function, |Vi?7p. This brings us to the re- 
quired saddle point x'^. After finding the location of the 
saddle point a;*, we calculate the exact direction of z ac- 
cording to Eq. (44), and the energy barrier according to 
Eq. (40). The direction z is then recorded and used as 
the initial guess for the next iteration in which the strain 
is further decreased. 
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FIG. 8: The probability for thermally induced plastic event as 
a function of 7 for different temperatures increasing from right 
to left. The system is of size A'^ = 484 and it has a mechanical 
instability at 7p = 0.0346. Inset: The mean distance 7p — (7) 
of the distributions for different temperatures. Here 7 = 2 x 
10^^ 
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FIG. 7: Scaling of energy barriers for various system sizes. 
The slope of the continuous lines is 3/2. 



In Fig. 7, we present typical AE{'-f) vs. 7p — 7 for 
systems of size N = 484,1024,2116,4096 and 10000. 
Clearly, as mentioned above, both the pre-factor and the 
range in strain for which the scaling law (41) holds are 
independent of N. 

Mechanical instabilities always involve the vanishing 
of the lowest eigenvalue of the Hessian; this observation, 
together with the above findings, imply that the process 
of loosing mechanical stability is reflected in the lowest 
eigenvalue of the Hessian only in some interval 7p — 7 
which is indeed system-size dependent. However, the 
same process is initialized at strain intervals that are sys- 
tem size independent, of the order of tenths of a percent 
in strain, see Fig. 7. This, in turn, implies that the strain- 
induced reduction of energy barriers is highly relevant for 
the discussion of thermally activated plasticity, for any 
system size. 



At higher temperatures this barrier can be overcome 
when 7 < 7p. To see the effect of temperature explicitly 
one introduces the probability to undergo a thermally 
activated plastic event at the strain value 7, denoted as 
Pa{-f,T) [14]. This probability was measured for a range 
of temperatures and is displayed in Fig. 8 for a system 
with N = 484 and an instability at 7p = 0.0346. This 
probability was measured at 7 = 2 x 10~^ by simulating 
a single elastic branch with randomized initial velocities 
corresponding to a given temperature. For each member 
of the ensemble we detect the value of the strain 7 at 
which the system leaves the basin of attraction of the 
athermal local minimum. For low temperatures the mean 
of the distribution is very close to 7p and the distribution 
is sharp. As temperature is increased the distributions 
move to the left, allowing a transition at lower values of 
7, with 7p — (7) being an increasing function of T. At 
about r = 0.1 for this model the distribution fiattens out, 
and for higher temperatures thermal noise is dominant 
over the mechanical instability. Using for this system 
cry « 0.5, /i « 15 and r* « 1 can estimate roughly 

AGtyp « 1. 



D. The third regime T > Tmax 

Having found an estimate for AGtyp at the present 
value of the strain rate, we can also estimate T'max(7) as 



AG. 



typ 



log((TY//i7r*) 



(45) 



Note that as a function of 7 this value changes apprecia- 
bly. For the present system this means that at tempera- 
tures higher than about 0.1 the plastic events are domi- 
nated by thermal noise. Since the glass transition in this 
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model is estimated to occur about about Tg « 0.46, we 
expect to have a sizeable region Tmax < T < Tg where a 
theory assuming that plastic events are uncorrelatccl and 
dominated by thermal noise might be a good model of 
the actual physics. 



V. CHOPPING OFF THE AVALANCHES: 
STRAIN RATE EFFECTS 

In this section we explain that increasing the strain 
rate results again (as for increasing the temperature) in 
turning the AQS anomalous stress fluctuation to a nor- 
mal process, but the physical mechanism is very different. 
Here the destruction of the correlated events is due to the 
forcing (by the faster strain rate) of simultaneous plastic 
drops, not letting them enough time to be correlated. To 
see this mechanism with clarity we need to expose a sec- 
ond length scale in addition to ^2 (T) that was defined in 
Eq. (34). This second length has to do with the strain 
rate 7. 



A. The typical length cissociated with strain rate 

We start by substituting Eq. (26) in Eq. (28) to obtain 
the scale s [15] . With A being the unit of length we write: 



cry A'' 



CMP 



(46) 



Consider next the rate at which work is being done at 
the system and balance it by the energy dissipation in 
the steady state, 



avW = (AC/)/rpi 



(47) 



where Tpi is the average time between plastic flow events. 
This time is estimated as the elastic rise time which is 



(Ao) 
M7 



ay 



(48) 



Next we note that Tpi decreases when A'' increases. On the 
other hand there exists another crucial time scale in the 
system, which is the elastic relaxation time Tei ~ L/c, 
where c is the speed of sound c = \f\jjp- Obviously 
this time scale increases with N like N^l'^ . There will be 
therefore a typical scale ^1 such that for a system of scale 
L = ^1 these times cross. At that size the system cannot 
equilibrate its elastic energy before another plastic event 
is triggered, and multiple avalanches must be occurring 
simultaneously in different parts of the system, each of 
which has a bounded magnitude. We estimate ^1 from 
Tel ~ Tpi, finding 



(a/c) 



UY 



A-^i 



A'^i 



(49) 



Using now the obvious fact that ~ (^/(A))'^ we 

compute 



6 
A 



GY>fi 



C 

A^ 



l/(l-/3d) 



(50) 



We observe the singularity for quasi-static strain when 
7 — > 0, where ^1 tends to infinity, in agreement with the 
results of quasi-static calculations. At low temperatures, 
before the thermal energy scale becomes important, the 
size of plastic flow events can be huge indeed. Note that 
there exists a difference between our estimate of ^1 and 
that of Ref. [19]. 



B. The effect of simultaneous plastic events: the 
Herschel-Bulkley law 

Due to the faster strain rates plastic events do not 
have time to cooperate and provide us with correlated 
avalanches. To exemplify the natural emergence of the 
length scale ^1 of Eq. (50) in this context, we mea- 
sure the difference between the mean flow stress mea- 
sured at the steady state (denoted (0^0)) and the con- 
figurational stress of the corresponding local minimum 
(denoted (uc)), see inset in Fig. (9). We expect this 
difference to go to zero in AQS conditions (since all the 
plastic events are discharged during the avalanches) and 
to increase with the strain rate. This difference is found 
by stopping an athermal simulation run at a given strain 
rate 7 (cf. Subsect. II B) and quenching the system to 
hit the local minimum where the configurational stress 
is measured. The quenching is achieved first by run- 
ning athermal molecular dynamics without increasing the 
strain for IOOt^ to allow for any plastic activity to com- 
plete, followed by a potential energy minimization. This 
procedure results in a stress drop as seen in the inset of 
the right panel of Fig. 9. The amount of stress drop 
is determined by the number of plastic events that occur 
starting at the point in time when the strain increase was 
stopped. 

Plotting individually ((Too) and (ctc) for one system size 
as a function of strain rate we note the relation of the 
present measurement to the Herschel-Bulkley law [21] 
which relates the mean flow stress to the strain rate, 



oy + 57 



(51) 



Here uy is the mean flow stress in the AQS limit and x is 
an exponent. This law is supposed to be N independent 
at large values of N . The power-law dependence implied 
by Eq. (51) is clearly seen in Fig. 9. 

We expect that the mean drop in stress between the 
flow stress and the configurational stress should be deter- 
mined by the ratio of the system size to the length scale 
^1 of Eq. (50). We thus propose a scaling function 



(52) 



10 




O ((Too) 



FIG. 9: The mean flow stress {a 00} and the mean configu- 
rational stress (ac) as a function of strain rate for a system 
with N = 1024. Inset: a typical drop in the strain when the 
system is minimized to the nearest local minimum. 
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FIG. 11: Left panel: the average difference between the steady 
state stress aoo and the corresponding stress at the local min- 
imum cTc for different system sizes as a function of 7 in three 
dimensions. System sizes are = 16000, (v), = 8000 (A) 
, N = 4000 iO) , N ^ 2000 (★), N = 1000 (0), N = 512 
(□). Right panel: Data collapse by re-plotting the same data 
as shown. The continuous line has a slope of unity. 
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FIG. 10: Left panel: the average difference between the steady 
state stress (Joa and the corresponding stress at the local min- 
imum (Tc for different system sizes as a function of 7 in two 
dimensions. System sizes are = 10000 (A) , = 4096 (O) 
, Af = 2116 (★), A^ = 1024 (0), N = 484 (□). Strain rates 
are 7 = 10~^ (full symbols), 7 — 10~* (dotted symbols) and 
7 — 10~^ (empty symbols). Right panel: Data collapse by 
re-plotting the same data as shown. The continuous line has 
a slope of unity. 

The dependence on iV* is called for by the clear N de- 
pendence seen in the left panel of Fig. 10. The exponent 
S and the asymptotics of the scaling function are deter- 
mined by (i) requiring the loss of the N dependence when 
7 — cx) and (ii) agreement with the Herschel-Bulkley law. 

We thus require gi{x) — > a;^ when x —>■ 00, and write 

S + C/d = 0, 1/(1- (3d) =x. (53) 

The best data collapse is obtained with C = 1; the cor- 
rectness of the scaling ansatz is shown in the right panel 
of Fig. 10 in which the data in the left panel is collapsed 
to a single function using the rescaling proposed by Eq. 
(52). The straight line in the right panel has a slope of 
unity to exemplify that C = 1- As a consequence of this 
result we find that 5 = —1/d (i.e. 1/2 in 2D and 1/3 in 
3D). Finally x = 1/(1 - (3d) which in 2D translates to 



FIG. 12: Left panel: The stress variance in athermal simula- 
tions (multiplied by A'^) for different strain rates as a function 
of the system size, see inset in the right panel. The change 
in scaling upon increasing the strain rate is obvious. Right 
panel: the data collapse in agreement with the scaling func- 
tion gi of Eq. (54) 

X = 3/7. 

These results continue to hold in three dimensions as 
well. We have measured /3 directly in 3D simulations and 
got the same number (3 « —2/3, cf. Fig. 2. Measuring 
the analog of Fig. 10 in 3D we obtained the data shown 
in Fig. 11. Indeed the result C = 1 continues to hold, 
consistent with S = —1/d and x in 3D being x = 1/3. 



C. Cross over in the fluctuation spectra due to 
strain rate 

As in the case of temperature, also strain rate destroys 
the highly correlated plastic events that are observed in 
AQS simulations. Again the relevant length scale should 
be £_i of Eq. (50), and we construct a scaling function 
to describe the change in the nature of the fluctuations 
using this same length scale. The data for the stress 
variance as a function of system size for different values 
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of the strain rate are shown in the left panel of Fig. 12. 
The change in scaling with increasing 7 is obvious. Here 
we propose a scaling function in the form N^^g{x), or 
explicitly 



(54) 



For 7^0 (equivalently ^1 ^ 00 or a; — >■ 0) we should 
recover the AQS results, requiring g{x) — >const. On the 
other hand when 7 is large, or a; — )■ 00, we should get 
uncorrelated fluctuations, {Sa"^) — )• N~^. This requires 
g{x) -J> a;^ when a; ^ 00 with C = 2{-20 - 1). These 
asymptotics, the cross over at a; = 1, and the excellent 
data collapse are all seen in the right panel of Fig. 12, 
adding full justification to estimate and the ramifications 
of the existence of the the typical scale ^1 . 

VI. THE UNIFYING SCALING THEORY 

In this section we provide a scaling theory that will 
unify the simulation results for all temperature and strain 
rates. Having two length scales at our disposal, we realize 
that the shortest of the two will be dominant at any given 
conditions (r, 7). We thus define ^ according to : 



T>T* 



(55) 



where T* is obtained by equating the two scales, i.e. 



~ e 



cry A''/ VA7 



(56) 



A very interesting and direct way of demonstrating the 
cross-over due to the combined thermal and finite strain- 
rate effects is provided by measurements of the variance 
of the stress fiuctuations as a function of the temperature, 
the strain rate and the system size. In Fig. 13 we dis- 
play 2D measurements of this quantity which is obtained 
by averaging the square of the microscopic stress fiuctu- 
ations in long stretches of elasto-plastic steady-states of 
the model described above at varying strain rates and 
temperatures as described in the figure legend. The vari- 
ance of the stress fiuctuations decreases as a function 
of A'', and in the left panel we multiplied the variance 
by N for better representation. Under quasi-static and 
athermal conditions the dependence is a power-law, cf. 
Eq. (29), with 6 ~ —0.4 in the present model. At higher 
temperatures and higher strain rate the data in Fig. 13 
indicate a clear cross-over to normal intensive stress fluc- 
tuations in which (Jcr^) ~ A^~^. 

As done before, to demonstrate unequivocally the 
cross-over that depends on both the temperature and the 
strain rate, we first need to separate the thermal from the 



mechanical contributions to (5(7^) cf. Eq. (32). Using the 
scale ^, we write the fiuctuations in the form 



V 



(57) 

The asymptotics of this function are again T{x) — > const 

for X — >■ and J^{x) — )■ a;^ for a; — )■ 00 as is born out in 
Fig. 13. 

VII. CONCLUDING REMARKS 

Wc have shown in this paper that the mechanical re- 
sponse of amorphous solids to external strains cannot 
be described with a theory that assumes a desert in 
which only thermal fiuctuations are important in trig- 
gering plastic events. Quite on the contrary, the problem 
abounds with rich and interesting correlated fiuctuations 
which are most spectacular in the AQS limit. This limit 
can be fully understood on the basis of the mechanical 
instabilities that are seen as eigenvalues of the Hessian 
matrix going to zero. There plastic events are highly 
correlated and appear in the form of system spanning 
avalanches. Increasing the temperature and strain rate 
has a great effect on these correlated plastic events, chop- 
ping them off and tending to make them normal in the 
limit of high temperature and/or high strain rate. We 
explained that the mechanisms for chopping off the cor- 
relations differ for temperature and strain rate, leading 
to a complicated cross-over behavior when T and 7 in- 
crease. There are three typical regimes, the AQS regime 
-^cross? intermediate regime ^c^oss *^ T <C ^max(O')) 
and a thermal regime when T > Tj^a^i'l)- The intermedi- 
ate regime is most interesting with its gradual change in 
the statistical physics of the system response and fiuctu- 
ations. Interestingly, we demonstrated that scaling con- 
cepts are of great importance in organizing the complex 
physics discovered and discussed above. All the effects of 
cross-over could be captured with the help of judiciously 
chose scaling functions whose existence means a high de- 
gree of predictability. By making measurements at some 
corner of the parameter space (T',7) we can predict the 
correct results for any other point in this parameter space 
with the help of the scaling functions presented above. 
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FIG. 13: Left panel: the measured stress fluctuations as a function of the system size for different temperatures and strain 
rates. Here the temperature ranges from T = to T = 0.15 and the strain rate ranges from 7 = 1 x 10~^ to 7 = 5 x 10~*. 
Right panel: data collapse as predicted by the scaling function Eq. (57) 
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